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Abstract 

A cosmological model in which primordial black holes (PBHs) are 
present in the cosmic fluid at some instant t = to is investigated. The 
time to is naturally identified with the end of the inflationary period. 
The PBHs are assumed to be nonrelativistic in the comoving fluid, to 
have the same mass, and may be subject to evaporation for t > to- Our 
present work is related to an earlier paper of Zimdahl and Pavon [Phys. 
Rev. D 58, 103506 (1998)], but in contradistinction to these authors 
we assume that the (negative) production rate of the PBHs is zero. 
This assumption appears to us to be more simple and more physical. 
Consequences of the formalism are worked out. In particular, the 
four-divergence of the entropy four-vector in combination with the 
second law in thermodynamics show in a clear way how the the case of 
PBH evaporation corresponds to a production of entropy. Accretion 
of radiation onto the black holes is neglected. We consider both a 
model where two different sub-fluids interact, and a model involving 
one single fluid only. In the latter case an effective bulk viscosity 
naturally appears in the formalism. 
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1 Introduction 



In studies of the early universe, the idea has sometimes been put forward that 
under these extreme physical conditions primordial black holes (PBHs) were 
present in the cosmic fluid during some brief time period. When considered 
in the comoving coordinate system, one may naturally take these PBHs to be 
nonrelativistic particles immersed in a cosmic fluid, most likely a radiation 
dominated fluid. The situation resembles that of a two-phase system of 
air bubbles in water, and similarly as in ordinary hydromechanics one may 
consider this physical system in one of the following two ways: either, as a 
system of two different fluid components in interaction, or, as an effective 
one-component fluid endowed with a bulk viscosity (under the assumption 
of complete isotropy the shear viscosity does not come into play) . 

The 1998 paper of Zimdahl and Pavon |l] (hereafter referred to as ZP; cf. 
also the more recent Ref. @) discussed a two- fluid model of the early universe, 
namely a nonrelativistic PBH component with particle masses lying within a 
narrow range, obeying the equation of state for dust, embedded in a radiation 
dominated cosmic fluid. The influence from the PBHs on the dynamics of the 
universe was described by the two-fluid model developed earlier by Zimdahl 
in Ref. §. 

A central ingredient of this theory was that the PBHs were allowed to 
be created, or annihilated, as long as the law of total energy conservation 
was not broken. One reason why this kind of model is physically appealing, 
is that the dynamics allows entropy to be produced. As is well known, one 
of the challenges of cosmological theory is to account for the large amount 
of specific entropy observed in the universe; the nondimensional entropy per 
baryon being about a ~ 10 9 ||. (It turns out that, if instead of a two-fluid 
model one chooses to work in terms of an effective one-component model, 
the interaction between the PBHs and the radiation fluid can be described 
by means of an effective bulk viscosity £; cf. Eq. (60) below.) 

The purpose of the present paper is to focus attention on the following 
point. ZP assumed that all the PBHs, with practically equal mass ttibh, 
were produced during a very short time interval prior to a a definite instant 
t = to- The magnitude of to was actually not given explicitly in Ref. but 
it seems reasonable to assume that to ~ 10 -33 s, i.e., inflationary times. (In 
an earlier work of Hayward and Pavon H the time to at which the PBHs 
were created was pushed further back in the history of the universe, to the 
order of a few Planck times.) For t > t , ZP assumed the evolution of the 
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PBH fluid to be governed by a fixed and negative production rate Tbh'- 

Tbh = --, (1) 

T 

r being the assumed constant life time of an evaporating PBH particle which 
at to has a given mass rriBH{to)- The history of the universe after t = to is 
thus in this picture determined by the Friedmann equations, together with 
the continuity of the PBH particles as determined by Eq. ([!]). The decay of 
the PBHs according to this picture, is analogous to the decay of an assembly 
of radioactive particles. Thus if one considers a comoving volume element, 
there is a definite probability for a certain fraction of the PBHs simply to 
disappear. 

It ought to be mentioned that the physical interpretation of Tbh is some- 
what problematic. It seems to be most natural to let this quantity refer to 
entire BH particles, so that Tbh means the specific fraction of BH particles 
disappearing, per unit time and unit comoving volume, from the cosmic fluid. 
"Specific" here means per particle. In general, if \I> denotes the production 
rate of particles per unit comoving volume, one has (nU^)-^ = where 

is the four- vector of the cosmic fluid and n the particle density (cf., for 
instance, Ref. |J). As — nT, one thus has the general equation 

(nU^ = nT, (2) 

showing explicitly how the quantity T refers to one particle. It may be 
instructive to compare with the interpretation of the Hoyle-Narlikar scalar 
C - field (cf, for instance, Ref. 0): in that continuous creation of 

matter is postulated to take place everywhere at any time; in the present 
case, a continuous destruction of BH matter may similarly be postulated to 
occur, during the life time r. 

One may however wonder: is this way of looking at the PBH evaporation 
physically correct? Imagine how the evaporation of the PBHs takes place: 
The black hole mass m,BH{t) changes with time, both because of absorption 
of electromagnetic energy (accretion), and because of Hawking radiation, but 
there appears to be no natural physical mechanism whereby some of the black 
hole particles should disappear suddenly. It may be argued that it would be 
more natural to assume that one black hole will not be able to disappear 
if not all the others do. And this is our main opinion to be forwarded in 
this paper: Assuming homogeneity and equal mass of all the black holes, we 
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argue that it is physically more plausible to assume that all the PBHs decay 
in exactly the same way, having the same life time r. Mathematically, this 
amounts to setting the negative production rate of the PBHs equal to zero: 

Tbh = 0. (3) 

The theme of our treatment in the following is thus to show how the universe 
develops after the instant t = to if the condition (3) is adopted, instead of the 
ZP condition ([]]). Therewith we avoid the problematic interpretation of Tbh 
altogether. No production or destruction of black holes is permitted. The 
formalism given by ZP will have to be changed and simplified. We will work 
out the changed formalism both for a two-fluid model, and for a one-fluid 
model. As already anticipated, the one-fluid model turns out to lead to an 
entropy production via an effective bulk viscosity, in a natural way. 

Before embarking upon an analysis of the two-fluid model, let us briefly 
consider the decay law for a single evaporating black hole. We assume that 
at time to all the PBHs with the same mass tube: (to) have been produced, 
and start evaporating. From Stefan-Boltzmann's law we have dmBn/dt = 
—aT BH A, where a is Stefan-Boltzmann's constant and A = l67rG 2 m BH is 
the surface area of a black hole with Schwarzschild radius 2Gm BB - Integra- 
tion with respect to mass from the instant £ onwards yields 

m BH (t) = (m 3 BH (t ) - 3C(t - t )) 1/3 , (4) 

where C is a constant. Since we associate all the PBHs with the same life 
time r, we require that m BB = at t = t + r. That determines the decay 
law; for later convenience we write it in the form 

m BH (t) -1 / t - t 
m BH {t) 3r V r 

The following point should be emphasized: we will ignore the effect of ac- 
cretion onto the black holes. This assumption puts a constraint on the mag- 
nitude of r, as can be seen from a simple argument: When a Schwarzschild 
black hole is immersed in an isotropic bath of photons with energy density 
p ra d far from the black hole, the accretion rate rh BH (t) is known to be f\ 

m BH (t) = 27tt G 2 m BH (t) p rad (t). (6) 
3 This argument ignores the universe expansion. 
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This follows from the capture cross section for photons from infinity being 
equal to cr phot = 27n G 2 m 2 BH (cf. Ref. [§). For the early times here in 
question, one has 

3H 2 1 
Pro, = g^, with H = -. (7) 

Now put t = 10~ 33 s. We obtain p rad = 1.9 x 10 54 GeV 4 , or 4.0 x 10 92 erg/cm 3 
which, when inserted into Eq. (6), is conveniently written as 

^ = 1.4 ( x 10 23 s-\ (8) 

mpi = G~ l l 2 being the Planck mass. This is thus the relative rate of 
mass increase due to accretion. Our condition about mass accretion being 
negligible, means that the expression (8) must be much smaller than the 
relative rate of energy loss due to evaporation. The latter is roughly equal 
to 1/(3t) according to Eq. (5). It follows that the condition on r becomes 

r « ( ^L) x 10- 24 s. (9) 

This condition will implicitly be assumed below. The greater the value of 
tubhi the stronger becomes the restriction on r. 

Let us also mention another point, related to our assumption that all 
black holes form almost simultaneously just prior to the instant t = t . As 
emphasized by Barrow et al. 0, in practice a hole of mass tubh cannot 
form until the cosmological horizon mass MnoRif) exceeds vtibh- Assuming 
radiation dominance, the horizon mass at time t is given by 

M HOR (t) = ^. (10) 

Thus, if we put t = 10~ 33 s, we see that the formed black holes must obey 
the restriction 

rriBH < 4 x 10 5 g. (11) 



2 The two-fluid system of PBHs and radia- 
tion 

In our picture there are two fluids present after the production time t : one 
PBH fluid whose decay in the comoving system of coordinates is given by 
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Eq. (5), and one radiation fluid whose description will be taken to be the same 
as given in Ref. ||]. Let us first consider the temperature behaviour of the 
PBH component. To this end we need a bit of relativistic fluid mechanics. We 
assume that the four- velocity of the cosmic fluid, U^ 1 = (U°, U l ), is common 
for the two sub-fluids. In a local inertial rest inertial frame thus U° = — Uq = 
1, U l = 0. The general expression for the fluid energy-momentum tensor is, 
for each sub-fluid component, 

T» v = pUW + (p + IL)hT, (12) 

where h^ u = + U^U U is the projection tensor, and II the extra pressure 
brought about by viscosity, by matter production, or eventually by both these 
effects together. Taking the four- divergence of Eq. (12) and multiplying with 
Uy, we have 

p + {p+ P + u)e = -u^. v , (13) 

where the scalar expansion is 9 = U^.^ = 3d/ a, a being the scale factor in 
the FRW line element 

ds 2 = _ dt 2 + fl2(t) ( + r 2 d tf\ (M) 

\ 1 — kr 2 J 

For the PBH sub-fluid, Hbh = 0. Neglecting also its pressure, pbh = 0, we 
get for the PBH component 

Pbh + PbhO = -UyT» v H . v . (15) 

From Eq. (2) we have in the local comoving frame 

ft + n6 = nT, (16) 

which for the PBH component implies 

n B H + n BH 9 = 0, (17) 

in view of our assumption (3). Putting psu = ^bh^bh we obtain from 
Eqs. (15) and (17) 

rtBHfrtBH = -UhTb H ;v- ( 18 ) 

4 As a digression, we note that according to Eq. (16) there could in principle be a short 
time period where T is positive, corresponding to a production of PBH particles. This 
production era must obviously be prior to tg. We do not model this era here. 
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Next, we consider the thermodynamical aspects of the system, starting from 
the thermodynamic identity (k B = 1) which holds for each sub-fluid compo- 
nent, 



nT a = p 



P + P . 
n, 

n 



(19) 



where a is the nondimensional entropy per particle. We regard n and T 
as independent variables, so that p = p(n,T). We differentiate this with 
respect to t, insert p from Eq. (13) and n from Eq. (16), and observe the 
useful relation 

p + p-Jg) = t(£) , (20) 



dn 



T 



dT 



which follows from the property of da being an exact differential [Q. We get 



T 
f 



dp) n T{dp/dT\ 



Tl6 + 



dp 
dn 



nT + 



This is the same form as given in Ref . , with the addition of the last term, 
which in turn reflects the fact that we are dealing with a non-closed system. 
When applied to the PBH component Eq. (21) yields, when T BH = 0, 



T, 



Pbh 



BH 



T, 



BH 



TT rpfiV 

ft BH-ui 



(22) 



where we have used that m BH = 1/(8ttGT bh ). From Eqs. (22) and (18) 



T, 



BH 



T, 



BH 



rriBH 
m B H' 



(23) 



This is actually the same as obtained in Ref. M. The relationship is inde- 
pendent of whether we put T B h = _ l/ r or Tbh — 0. It might perhaps seem 
surprising at first sight that the temperature in the black-hole fluid is not 
influenced by an eventual (negative) production rate of black holes. How- 
ever, the present model is based upon a perfect fluid assumption; therefore 
one should not expect that an existing black hole "feels" if another black 
hole is created or annihilated. The temperature Tbh is associated with each 
black hole separately, and cannot be thought of as a conventional common 
temperature for the whole fluid. 

Consider now the entropy four- vector, = noU^ . This expression holds 
for each of the components separately. Taking the four- divergence of it 
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we obtain, when using Eqs. (13), (19) as well as the production equation 
= nV, 

S^. fl = ~e-^nT-^U fl T^, v , (24) 

with p = (p + p)/n — To being the Gibbs potential per particle. As S' M ;(U = 
nib + aT) in the local rest frame, we can alternatively express Eq. (24) in 
the form 

^-^"-l^ r -^v < 25 » 

Again, this is an appropriate generalization of the corresponding expression 
given in ||. When applied to the PBH component Eqs. (24) and (25) yield, 
in view of Hbh = 0, T BH = as well as Eq. (18), 

Sbh-.h = n BH^ — , (26) 

or 

(Tbh = 7f — ■ I 2 ') 

J-BH 

Here Eq. (27) has the form that we would expect. 

We move on to the radiation sub-fluid. With H ra d = 0, and since generally 
l^rad = 0, we first obtain when using Eq. (24), 

CM _ _ ^ 77 rpflV _ -I jj rpiLV ,cyQ\ 

°rad;iJ, ~ rp ° M J rad;u ~ rp U M J BH;V l Z ° ) 

J- rad J- rad 

where in the second equation we have exploited the zero-divergence condition 
for the total energy-momentum tensor T^ d + T BH -. Thus 

SU, = ~ n BH ^ (29) 

rad 

in view of Eq. (18). From the specific entropy equation (25) we obtain, when 
setting p ra d = \p T ad and assuming entropy conservation, b ra( i = 0, 

3 pbh rriBH , Qn x 

L rad — — ■ 1<JUJ 

4 Prad m BH 

Thus, whereas we could above put Tbh = in the case of pure black- 
hole evaporation, this is not the case for Y ra d. When black holes evaporate 
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if^BH < 0), one gets a positive value of T rad and radiation particles become 
produced, as one would expect. 

Taking the four- divergence of the total entropy four-vector we get from 
Eqs. (26) and (29) 

5 % = S BH;n + S rad^ = PBH ~ — [ Tf, Ti ' ( 31 ) 

m BH \J-BH J- rad J 

This is an important result. The second law of thermodynamics requires S 7 ^ 
to be non-negative, so that the case of pure BH evaporation we must have 
Tbh > Trad- Black hole evaporation is connected with entropy production. 
If the converse were true in some time period, Tbh < T ra d, then the second 
law would require that ttibh > 0. The black holes would accordingly increase 
in mass during such a period. 

For comparison we write down some expressions obtained on the basis of 
assuming instead T B h = _ l/ r > as in the ZP paper. Instead of Eqs. (30) and 
(31) one gets 

r 3 pen ( 1 rn BH \ , , 

i rad\ZP - T 1 , (32J 

4 prad V r m BH J 



s^ ZP = -esa(^L') +PBH 'HsiL('-^). (33) 

r \m B H T rad J men \T B h T rad J 

The structure is here seen to be more complicated, especially for the entropy 
four-divergence. 

We will need also the expression for the relative rate of change of the 
temperature for photons. From the general expression 

^ = -\{e-Trad) (34) 
■J-rad <J 

we get, assuming Tbh = 0, 

Trad = J^( d+ 3p m mBH\ ^ 



Trad 3 \ 4 Prad mBH ) 

whereas assuming Tbh = _ l/ r we S e ^ 

T ra d\ZP 1 



■rad\ZP 



e+ 3f>BH_ fl m BH \ 

4 prad V T ITIbh) 



(36) 
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2.1 What happens for t>to? 



We are now able to give a physical interpretation of the two-component cos- 
mic fluid near the production time t = t . As mentioned above, we associate 
to with the end of the inflationary era. At this instant we will moreover make 
the explicit assumption that the sub-fluids have the same temperature. That 
is, the PBHs are taken to be created with the same temperature as the radia- 
tion field which they were created from. Since the two fluid temperatures co- 
incide, Eq. (31) tells that no entropy will be produced at t . Recall that at t , 
the second law of thermodynamics permits black-hole mass accretion as well 
as evaporation to occur. The first case corresponds to tHbh > 0, Tbh < T rad , 
the second case corresponds to rh BH < 0, T BH > T rad . 

At this point we have to make a digression, to check that our assumption 
about equality of the two sub-fluid temperatures at t — to makes physi- 
cal sense. Evidently, the produced black holes have to be smaller than the 
universe itself. At time t the particle horizon R p (t), i.e., the radius of the 
observable part of the universe, is 



where we have used that a(t) oc t 1/2 . At t = 1(T 33 s thus R p = 6 x 1(T 23 
cm. As mentioned above, p ra( i = 4.0 x 10 92 erg/cm 3 , and the temperature 
T rad can be found from the equation p rad = a rad T^ ad , where a rad = 7.56 x 
10~ 15 erg/ (cm 3 K 4 ) is the radiation constant. We find T rad = 4.8 x 10 26 K, or 
4.1 x 10 13 GeV. The Schwarzschild radius of a black hole is tbh = 2Gm B H = 
1/(AttTbh) and we get, when assuming T rad = T B h, f bh = 3.8 x 10~ 29 cm. 
This is seen to be much smaller than R p . There is thus ample room for the BH 
gas in the universe, showing that our assumption about equal-temperature 
of the sub-fluids at t — t is consistent. 

Assume next that rfiBH = at t . Then no process at all, except from 
the expansion of the universe, will run. According to Eq. (23), T B h will then 
stay constant, whereas T rad will drop due to the expansion. That means, for 
t > t one has T rad < Tbh, and an evaporation process has to start, for the 
second law of thermodynamics not to be violated. 

If on the other hand an evaporation process is already going on at to, the 
black hole temperature T B h will have to rise according to Eq. (23). Equation 
(35) then tells that the radiation temperature T rad may rise or fall depending 
on the ratio Pbh/ Prad- If the process of pure evaporation is to be maintained, 




(37) 
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the radiation temperature must not rise faster than the black-hole tempera- 
ture. From Eqs. (23), (35), and (5), this leads to the following condition at 
to: 

/3<4(fl r + l), (38) 

where we have defined (3 by 

P = P ~^y 09) 

PradV'O) 

9q meaning the scalar expansion at to- If the inequality (38) is broken, the 
radiation temperature will at some time exceed that of the black holes, and 
the evaporation process will no longer be thermodynamically allowed. The 
maximum value of (3 represents the largest density ratio for which evaporation 
can take place. The possible values of (3 depend on the expansion and the 
life time of the black holes. The larger the value of r, the larger may the 
initial abundance of the black holes be. This is physically reasonable: a long 
life time means that the black hole is large and evaporates less intensively 
than a small black hole. We should therefore expect that large black holes 
contribute less to the heating of the radiation sub-fluid than small black holes 
do. Recall, however, the restriction on the magnitude of r that lies in Eq. (9). 

It is instructive to get an idea about the typical magnitude of the quantity 
9 t. Let us assume, as an example, that black holes of mass m BH = 100 g 
0.5 x 10 7 mpi are formed at t — 10~ 33 s. (This mass satisfies the condition 
(11).) From Eq. (9) we then have r <2x 10~ 31 s, in order to permit neglect 
of accretion. The scalar expansion is 9 = 3H = 3/(2£) = 1.5 x 10 33 s _1 , so 
that we obtain 9 t <C 300. 



2.2 Reheating temperature 

We assume henceforth the process of pure evaporation, so that for t > t the 
black-hole temperature increases and stays above the radiation temperature. 
The decay rate, Eq. (5), will then hold for the entire period t > t . 

As the second law of thermodynamics limits the reheating, a reheating 
temperature may be defined. It is the maximum temperature obtained 
for the radiation sub-fluid, and corresponds to T rad = 0. Using the rate 
formula (35), the equation of state p rad = Sn rad T rad , and Eq. (5), we find , 
when assuming the reheating to take place immediately after t = t : 

T reh = J_ n BH{to) m BH (t ) 

rad 12 n rad (t ) 9 t ■ 1 ) 
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By contrast, ZP obtained a value four times greater: 

rpreh a rpreh ( A~l\ 

1 rad\ZP — q:1 rad- K^ 1 ) 

This difference is physically understandable from a comparison between 
Eqs. (32) and (30). In the ZP case there is one extra term proportional 
to 1/t, contributing to the reheating. One would therefore expect that a 
higher amount of energy exists in the radiation field before the second law 
of thermodynamics forces the heating of the radiation to terminate. 

2.3 Equilibrium solutions 

Some years ago Barrow et al. |9| found cosmological solutions for which the 
ratio of the energy density of the black hole component to the energy density 
of the radiation component remained equal to a constant n, 

PBH = UPrad- (42) 

Solutions of this kind are called equilibrium solutions. 

The question naturally emerges: are such solutions possible in our model? 
To investigate this point, let us first observe that Eqs. (15) and (18) provide 
us with the following equation in the case of pure evaporation: 

Pbh + PbhO = n BH m BH . (43) 

The corresponding equation for the radiation component reads 

4 

Prad + ^Prad 9 = n BH m BH . (44) 

Inserting Eq. (42) into Eq. (43) and making use of Eq. (30) we obtain 

r md = -~ K + — V ( 45 ) 

4 V Prad J 

With p ra d = 3n ra d T ra d the last term to the right becomes 

Prad ^rad . Trad ^, „. , . 

= + TfT- = -{Trad-0), (46) 

Prad T^rad -1 rad " 

in view of Eqs. (16) and (34). Altogether we thus obtain, when Ybh = 0, 

Trad = \^-e. (47) 
4 K + 1 
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This shows that for an equilibrium solution to exist, T rac i has to be pro- 
portional to 9. Since T rad = —{?>K/4)rn BH /m BHl this implies in turn that 
ttibh I m BH has to be proportional to 9, or 

eocfl-t^-Y 1 , (48) 



according to Eq. (5). A relationship of this form is physically unlikely. Such 
a time dependence for 9 is quite different from that given, for instance, in 
|T0| (for a viscous universe). We conclude that equilibrium solutions do not 
lead to a realistic space-time dynamics when V bh — is assumed. The 
reason for this deviation from the theory of Barrow et al. is easy to trace 
out: these authors consider a very different kind of PBH model, in which the 
initial number of holes between m and m + dm is taken to be proportional to 
m~ n dm, with n (the spectral index) lying between 2 and 3. This evidently 
contrasts our adopted model where all the PBHs have initially the same 
mass m BH (to). We have no reason to expect that the consequences of two so 
different models are in agreement. 

The possibility of equilibrium solutions was investigated also by ZP: they 
found that by assuming Ybh = — 1/t, equilibrium solutions were mathemat- 
ically impossible. The qualitative agreement between the consquences of the 
ZP theory and the closely related present theory is of course what we would 
expect. 

2.4 Time dependence of the densities 

Integrating Eq. (43) from t to t, inserting 9 = 3d(t)/a(t) and using Eq. (5), 
we get 

PBH (t) = PBH (t ) (7) 3 ( 1 "^) l/3 - ( 49 ) 

Equation (44) for the radiation now becomes 

a pBH{to) fa \ 3 ( t-t V 2/3 fnn\ 

Prad + 4- prad = ~ I — J U • (50) 

a 3t V a J \ r J 



We multiply this equation with a 4 and integrate 

Pbh{Po) a o [ l _ fJ j- s ( \ t' — t \ 2 ^ 3 
a 4 3r a 



Prad{t) = prad{to)— H J / a(t ) [ 1 — — ) dt . (51) 



to 



13 



Knowledge about p rad (t ) and Pbh{Po) plus solution of the field equations 
for a(t) thus determine the dynamics of the entire system. In particular, for 
a flat FRW universe the first Friedmann equation yields 

a 2 8nG , . 

— = 7, \PBH + Prod)- (52) 
CL O 

Equations (49) and (51) are simpler than the corresponding equations ob- 
tained by ZP by assuming nonvanishing Tbh- 



3 On an effective one-component model 

Instead of working in terms of a two-component fluid model let us now try 
to formulate an effective one-component model. As mentioned earlier, this is 
analogous to treating water with air bubbles as one single fluid. The trans- 
formation from a two-fluid description to an effective one-fluid description, 
in the form worked out by Zimdahl ||, is not directly applicable here since 
there exists a temperature difference between the two components: the evap- 
orating black holes become hotter, while the radiation may cool down due to 
the expansion of the universe. We may nevertheless introduce an "equilib- 
rium" temperature T, but evidently so only in a formal sense except in the 
initial stages near t — to- 

As we will see, the lack of a physical equilibrium temperature does not 
prevent us from constructing in some sense an effective one-component the- 
ory. A key element in the formalism is that the evaporation process is repre- 
sented in terms of a viscous pressure. As above, we assume henceforth that 
Tbh = 0. 

Following the ZP paper we start from the Gibbs equation for a two-fluid 
model close to equilibrium: 

Tda = d (J^j + pd f -J - (u BH - p rad )d (^^j , (53) 

where T is the equilibrium temperature for the system of a conventional 
fluid and a PBH component. In our case the equations of state are p = 
Prad(n rad ,T) and p = p(n,n BH ,T), with n = n BH + n rad . We define the 
"equilibrium" temperature T through the equation 

PBH{n B H, T B h) + Prad{n rad , T rad ) = p(n, n B n, T). (54) 
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Since moreover [i rad = we have 



n/i = n BH fi BH , nT = n rad T rad , (55) 

na = n BH &bh + n rad a rad , (56) 

whereas upon use of pbh = 0, a BB = ^Gm BH , T B u = l/(SnGm B H) in the 
relation p BH = m BH - a BH T BH we get 

HBH = -jTYlBH, W = t^Pbh- (57) 

The energy-momentum tensor is given by Eq. (12) as before, whereas Eq. (13) 
because of conservation of total energy is replaced by 

p+(p + p + Tl)6 = (58) 

in the effective one-component description. From Eq. (53) we now obtain, 
after a brief calculation observing that {n/n rad )Y = T rad , 



na 



.«n + ( BM ,^-2±£)r. (59) 



By means of Eqs. (57) we get the following expression for the four-divergence 
S^-,n = n(& + oT) of the effective entropy four-vector: 

S^^—U. (60) 

This formula has a remarkably simple appearance, since the effective pro- 
duction rate T is absent. The whole of the entropy production is attributable 
to the excess pressure II, which for a conventional viscous fluid is expressed 
as II = — (9, ( being the bulk viscosity. Thus as far as entropy production 
is concerned, in our model the universe, composed of radiation and non- 
relativistic evaporating PBHs at a higher temperature than the radiation, 
can be described in terms of an effective bulk viscosity. 

The form of Eq. (60) relies upon the condition T BH = 0. The correspond- 
ing Eq. (68) in the ZP paper is more complicated. 

We may also compare the two expressions obtained for the production 
of entropy. We recall that in the two-component model, expression (31) was 
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derived on basis of the assumption about black hole evaporation. Comparing 
this expression with expression (60), we get 



where Eq. (30) has been taken into account. As we noted earlier, pure evap- 
oration corresponds to T BH > T rad , so that II has to be negative according 
to Eq. (61) since T ra( i > 0. (In other words, ( = — U/6 has to be positive, as 
it should.) 

Summing up, we see that the legitimacy of the one-component model 
rests upon the possibility that we have to put the temperatures of the two 
sub- fluids equal, at t = t . Although some of results of this theory are 
simple, in particular the expression (60) for the entropy, we think on the 
whole that the one-fluid model is from a physical point of view inferior to 
the two- fluid model considered earlier. The "equilibrium" temperature T 
defined in Eq. (54) is after all, similarly as in the ZP paper, introduced in 
a formal sense, so that the main physical relevance of the one-fluid model 
appears to be in the period just after t = t , when the two temperatures have 
not got the opportunity to separate significantly. 

4 Concluding remarks 

Our main new contribution in the present paper has been to provide an anal- 
ysis of cosmological BH theory when the production rate Tbh is assumed to 
be zero throughout. This assumption is based upon both physical and math- 
ematical considerations. Physically, the interpretation of a production rate 
in the form T BH = — 1/r, as assumed in the ZP paper M, is problematic (we 
made an attempt to formulate a natural interpretation in Sec. 1), whereas 
mathematically the condition Tbh = leads to considerable simplifications. 
In particular, the expression (31) for the four-divergence of the entropy four- 
vector fits in with the second law of thermodynamics in a simple way, showing 
how the evaporation of the PBHs leads to a positive entropy production. In 
principle, a process of this sort may even be one of the reasons for the re- 
markably high specific entropy, a ~ 10 9 , in the universe. 

Another basic assumption in our theory is that all the PBHs have the 
same mass. Physically, this is consistent with assuming that all the PBHs 
were formed during a short time interval just preceding the instant t = t Q , 




(61) 
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in the early universe. The time to is naturally taken to be at the end of the 
inflationary era. Immediately before to, the PBHs are assumed to be formed 
at the radiation temperature. Immediately after to, the PBHs may evaporate 
freely and the thermal equilibrium between black holes and radiation is lost. 

When the PBHs are evaporating for t > t , the radiation temperature 
T ra( i cannot rise faster than then black. hole temperature The reheating tem- 
perature is given by Eq. (40), and is four times higher than the reheating 
temperature (41) obtained in the ZP formulation. 

It should be noted that our theory neglects the accretion effect for the 
black holes. The mathematical condition for this is given by Eq. (9). 

When considering instead an effective one-component model, the simple 
appearance of the four-divergence formula (60), again being a consequence of 
Ybh = 0, should be noticed. This formula shows how the entropy production 
can be ascribed to the existence of an effective bulk viscosity ( = — 11/9 
in the mixed, common fluid. However, on the whole, the one-component 
model seems to be less motivated physically than a two-component model. 
The difference in temperature between the BH component and the radiation 
component signalizes a physical instability which is best handled physically 
when one regards the cosmic fluid as two separate sub-fluids. 

Finally, to widen the scope of applicability of the above formalism,we 
mention another interpretation of it that also appears to be physically pos- 
sible. Instead of requiring T BH to be strictly zero, as we have done above, 
one can assume that this requirement is only as approximation to the real 
physics. In that case, the equation Ybh = would correspond to a neglect 
of possible remnants of the PBHs. This would restrict the above formalism 
to hold only in the restricted time range t — t < r. 
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